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We describe a quasiperiodic optical lattice, created by a physical realization of the abstract cut- 
and-project construction underlying all quasicrystals. The resulting potential is a generalization 
of the Fibonacci tiling. Calculation of the energies and wave functions of ultracold atoms loaded 
into such a lattice demonstrate a multifractal energy spectrum, a singular continuous momentum- 
space structure, and the existence of controllable edge states. These results open the door to cold 
atom quantum simulation experiments in tunable or dynamic quasicrystalline potentials, including 
topological pumping of edge states and phasonic spectroscopy. 
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I. INTRODUCTION 

Quasiperiodicity has a profound impact on electronic 
structure, playing a role in phenomena ranging from the 
quantum Hall effect to quasicrystalline ordering. How¬ 
ever, the formation, stability, excitation, and electronic 
structure of quasiperiodically ordered systems remain in¬ 
completely understood. Open questions include the na¬ 
ture of electronic conductivity or diffusivity, the spec¬ 
tral statistics, the nature of strongly correlated magnetic 
states on a quasicrystalline lattice, topological proper¬ 
ties of quasicrystals, and even the shape of the electronic 
wave functions [iHH]. 

The exquisite controllability of cold atoms makes them 
a natural choice for experimental investigation of the 
open questions regarding quasiperiodicity. Unique fea¬ 
tures of such experiments would include precisely variable 
quasiperiodic parameters, tunable interactions, bosonic 
or fermionic quantum statistics, and the ability to study 
dynamical phenomena (in modulated or quenched sys¬ 
tems, e.g.). Numerous theoretical proposals have ex¬ 
plored the rich physics of quasiperiodically trapped cold 
atoms [EHZI]. However, with the exception of some 
early experiments on non-degenerate atomic gases in 2D 
quasiperiodic lattices [2H1 [53] j the dominant application 
of quasiperiodic or incommensurate potentials in cold 
atomic physics thus far has been as a convenient proxy 
for disorder [301 ini e.g.]. The realization of tunable 
quasicrystalline potentials for cold atoms would open up 
a broad range of exciting experiments, complementary 
to those possible with synthesis and characterization of 
solid-state or photonic quasicrystals. 

In this paper, we describe and elucidate the proper¬ 
ties of a “generalized Fibonacci” optical lattice which 
creates a dynamically tunable family of ID quasiperi¬ 
odic structures. This lattice physically realizes the ab¬ 
stract cut-and-project construction which underlies all 
quasicrystals. Every quasiperiodic tiling can be defined 
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FIG. 1. (Color online) The generalized Fibonacci optical lat¬ 
tice. Top: Diagram of cut-and-project construction of the 
generalized Fibonacci lattice. A 2D strip at a particular angle 
a is projected down to a line. If tan(a) = l/r = 2/(l-|- \/5), 
this results in the Fibonacci tiling itself; a different irrational 
slope creates a different ID quasicrystal. Bottom: Calcu¬ 
lated potential of a Fibonacci optical lattice with parameters 
discussed in text (red is deeper), showing Fibonacci sequence 
rlrrlrlrrlTTlrlr... of lattice spacings. 


as a projection of a cut through a lattice which is fully 
periodic but exists in a higher dimensional space [H [35] . 
For example, the Fibonacci tiling is a projected ID cut 
from a 2D square lattice, and the Penrose tiling can be 
constructed as a projected 2D cut through a 5D peri¬ 
odic lattice. In these and all other quasicrystals, the hid¬ 
den degrees of freedom in the higher-dimensional space 
can give rise to phenomena such as phasonic excita¬ 
tions and Bragg diffraction with unconventional symme¬ 
try [5] 1331 133 . Generalized Fibonacci optical lattices 
will provide a flexible platform for realization of tunable 
quantum quasicrystals, and should enable direct exper¬ 
imental investigation of questions inaccessible to exper¬ 
iments on static, non-tunable, non-interacting systems. 
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Specific topics of interest include studies of edge states, 
adiabatic quantum pumping, multifractal energy spectra, 
phasonic spectroscopy, dynamical signatures of many- 
body localization, and transport in quasicrystals. 


II. CONSTRUCTION OF GENERALIZED 
EIBONACCI LATTICES 

The generalized Fibonacci optical lattice is constructed 
as a direct real-space realization of the cut-and-project 
procedure, by intersecting an elongated optical trap at a 
tunable angle with a large-period square lattice, as dia¬ 
grammed in Fig. The resulting potential is the sum of 
two simple red-detuned trap potentials: the large-period 
square optical lattice, with potential 

VLix,y) = -Al sin^ (^) ’ 

and the elongated cutting beam at an angle a to the x 
axis, with potential 


Vc{x,y) = - Ac exp 


-2 


—x sin(a) -I- y cos(a) 


Wo 


Here Ac is the depth of the square lattice, Ai/2 = a 
is the lattice constant, Ac is the trap depth of the 
cutting beam, wq is the beam waist, and we have as¬ 
sumed a Rayleigh range long compared to the trapping 
region. The natural energy scale is the recoil energy 
Efi = h?/2m\\. The total potential is then U{x,y) = 
Vl{x, y) + Vc{x, y, a). This potential is shown in Fig. 
for tan(Q!) = 2/(1 -|- y/h), and Ac = lOA^. In order 
for the total physical potential to be a good approxima¬ 
tion to a true cut-and-project potential, Ac/Ac must 
be sufficiently large that we can spectrally distinguish 
states along the cutting beam from transversely-extended 
states. To preserve the ID character of the potential, the 
waist of the cutting beam should not be large compared 
to the lattice constant a; however, the calculations pre¬ 
sented below indicate that the potentials retain many of 
their interesting quasiperiodic properties even if this con¬ 
dition is violated. A perturbative treatment elucidating 
the quasiperiodic character of the effective potential ap¬ 
pears in the appendix. This trap construction can be 
generalized in a straightforward way to 2D quasiperi¬ 
odic traps, via intersection of a light sheet with a 3D 
large-period lattice. This direct experimental realization 
of the simplest quasiperiodic lattices is also intrinsically 
tunable: variation of the intersection angle a tunes the 
properties of the resulting potential, generating different 
members of this family of quasicrystals, and variation of 
the offset transverse to the cut beam axis drives phasonic 
degrees of freedom [HH] . 

We now briefly discuss the practical optics required to 
realize such a potential. If the cutting beam is produced 
by focusing a gaussian beam of initial diameter D with a 


lens of focal length F, the number of lattice sites in one 
Rayleigh range of the beam is given approximately by 


All 


16 Ac /FV 
TT Al 'vfJ ’ 


where Ac is the wavelength of the cutting beam and Ac 
is twice the lattice period. Even if these traps are pro¬ 
duced by the same laser, Ac/Ac can be varied by using an 
angled-beam lattice configuration. The number of lattice 
sites spanning the cutting beam width is 




8 ^ /FA 
^ Ac 'vfJ ’ 


so the aspect ratio of the full trap is A||/Aj_ = 2F/D. 
With typical values of F and D, one can then real¬ 
ize a range of generalized Fibonacci traps, with widths 
ranging from less than a lattice constant to many lat¬ 
tice constants. The ends of such a trap can be de¬ 
fined for example by tightly-focused blue-detuned light 
sheets. The intersection angle a is most easily tuned by 
rotating the lattice itself; for angled-beam lattices cre¬ 
ated using a diffractive optical element, this could be 
straightforwardly achieved with a single rotation stage. 
As with ordinary optical lattices, adiabatic loading of 
cold atoms into a Fibonacci-type lattice would be ac¬ 
complished starting from the elongated optical trap by a 
slow turn-on of the lattice potential. 

Using this trap geometry, one can construct a contin¬ 
uous family of quasiperiodic tilings of the line by placing 
the cutting beam at any irrational slope. If the angle 
of intersection a satishes the relationship tan(a) = 1/r 
where r is the golden mean (1 -|- -\/5)/2, then the re¬ 
sulting potential will approximate the Fibonacci tiling; 
this is the reason we refer to the family of potentials as 
“generalized Fibonacci” lattices. This one-dimensional 
structure tiles the line quasiperiodically, exhibits sharp 
diffraction peaks, and can also be generated algebraically 
using the substitution rule r — >■ rl, 1 —r, which gives 
rise to the sequence (I, r, rl, tIt, rlrrl, rlrrlrlr, 
tIttItIttIttI...). As the bottom panel of Fig. 
demonstrates, the energy minima of the Fibonacci op¬ 
tical lattice are spaced according to the Fibonacci tiling. 
Because of the inflation symmetry associated with the Fi¬ 
bonacci tiling, if the width of the cut-out strip is reduced, 
then the resulting one-dimensional projection will simply 
be an expanded and displaced version of the Fibonacci 
tiling [36] . In the generalized Fibonacci optical lattice, as 
the width of the Gaussian cutting beam is increased, the 
potential no longer approximates a one-dimensional pro¬ 
jection, but remains quasiperiodic. The connection to 
a mathematically exact cut-and-project lattice and the 
quasiperiodicity of the potential emerge clearly from a 
perturbative treatment, which also makes plain the con¬ 
nection to closely related systems such as the Aubry- 
Andre model. Details of such a perturbative treatment 
appear in the appendix. This general optical technique 
for construction of a family of quasiperiodic lattices and 
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FIG. 2. (Color online) Energies and wave Inactions in a tun¬ 
able Fibonacci-type potential, a: A portion of the energy 
spectrum of the generalized Fibonacci optical lattice as a func¬ 
tion of cut angle. Here a = 1, Ac = 5Ac, wq = 1. 

Color of points corresponds to center-of-mass of probability 
density, to enable identihcation of edge states, b: 2D probabil¬ 
ity density versus position at the indicated cut angle and en¬ 
ergy (a typical bulk state). Red regions correspond to higher 
probability density, c: 2D probability density versus position 
at the indicated cut angle and energy (a typical edge state). 
Both b and c show a region of 4 by 200 lattice constants. 



FIG. 3. (Color online) Tunable quasicrystals in momentum 
space. Main lower panel shows Fourier transform of the 
ground state probability density along the direction of the 
cutting beam, as a function of cut angle a. Note the logarith¬ 
mic scale on the colorbar. Upper panels show log of Fourier 
transform amplitude versus wave vector at a cut angle of 7r/4 
(top) and tan~^ (2/(l-|--\/5)) (bottom), with identical axis lim¬ 
its. Dashed lines show expected wave vectors of Fourier peaks 
of the potential based on the second-order perturbation the¬ 
ory described in the appendix. Calculations were performed 
on a mesh 4 by 200 lattice constants in size. Faint vertical 
lines in main plot are hnite-size edge effects. 


their rational approximants is the first main result of this 
work. 


III. WAVE FUNCTIONS AND ENERGIES IN 
GENERALIZED FIBONACCI LATTICES 

The second main result of this work is the calcula¬ 
tion of the energy spectra and wave functions of non¬ 
interacting atoms trapped in this family of tunable qua¬ 
sicrystalline potentials. These calculations demonstrate 
the utility of generalized Fibonacci optical lattices as 
a tool for the investigation of quasiperiodic quantum 
phenomena. To determine the energy spectrum of the 
physically-realized trap, we solve the two-dimensional 
single-particle Schrodinger equation on a mesh with spac¬ 
ing much smaller than a lattice constant. This approach 
avoids simplifications inherent in the tight-binding ap¬ 
proximation, and makes closer contact with experimen¬ 
tally realizable traps. We do not use periodic boundary 
conditions, both for more direct comparison with real ex¬ 
periments and so as to accurately model the existence of 
edge states at the ends of the quasicrystal. Energy eigen¬ 
values as a function of cutting beam angle are shown in 
Fig.i The calculated spectrum has a complex multifrac¬ 
tal appearance m- Notable features include a hierarchy 
of minigaps which disperse as the angle is varied, a non¬ 
accidental resemblance to the Hofstadter butterfly, and 


the existence of isolated states in the gaps. We find that 
the qualitative structure of the energy spectrum remains 
the same if the waist of the Gaussian beam is increased to 
several times the size of the lattice constant. The resem¬ 
blance to the Hofstadter butterfly is to be expected, given 
the recent demonstration that the generalized Fibonacci 
quasicrystal and the Harper model of high-magnetic-field 
2D integer quantum Hall states are topologically equiv¬ 
alent [551 [55] . The intersection angle a of the Fibonacci 
quasicrystal plays a role analogous to that of the modula¬ 
tion period of the Harper lattice, or the effective magnetic 
field in the quantum Hall system. 

The wave functions of atoms in generalized Fibonacci 
optical lattices also possess unique characteristics. As 
the cut angle a is varied, the Fourier transform of the 
projected spatial probability density of the ground state, 
plotted in Fig. shows for irrational tan(a) a rich sin¬ 
gular continuous structure characteristic of quasiperiodic 
structures. This property is the ID analogue of the for¬ 
bidden Bragg diffraction patterns by which 3D quasicrys¬ 
tals were first discovered m, and recalls the definition 
of a quasicrystal as a structure which produces a sharply 
peaked diffraction pattern but lacks translational symme¬ 
try. The singular continuous nature of the spectrum also 
emerges naturally from a perturbative treatment of the 
effective potential (see appendix for details). Rational 
tan(a) = p/q produces a crystalline superlattice with a 
periodicity which depends upon q. The real-space struc¬ 
ture of a typical squared wave function in a Fibonacci 
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optical lattice is shown in Fig. In addition to ex¬ 
tended bulk states, isolated states traversing the band 
gaps are visible in Fig. [§L. Fig.[^ shows the squared 
wave function of a typical gap-traversing state, located 
inside the lowest energy band gap, and demonstrates that 
the wave function is localized towards one edge of the lat¬ 
tice. These states, which can occur at both rational and 
irrational cut slope, are interesting candidates for realiz¬ 
ing topological pumping iHiniiia- 


IV. SELECTED APPLICATIONS: 
TOPOLOGICAL PUMPING AND PHASON 
SPECTROSCOPY 


Topological pumping is possible because the wave func¬ 
tions and energy spectra of the generalized Fibonacci op¬ 
tical lattice depend in a non-trivial way on the offset of 
the cutting beam with respect to the lattice. In the ter¬ 
minology of quasicrystals, this offset is a phasonic degree 
of freedom. Just as phonon modes arise from discretely 
broken real-space translation symmetry, phason modes 
arise from broken translation symmetry in the higher¬ 
dimensional space from which a quasiperiodic lattice is 
projected [33l |35l |4T] . In a Fibonacci-type optical lattice, 
this corresponds to symmetry under relative translation 
of the cut beam and the lattice in a direction transverse 
to the cut beam. A visualization of the effects of contin¬ 
uous adiabatic phasonic driving in the Fibonacci optical 
lattice is shown in Fig. As the offset of the cutting 
beam is varied from 0 to I lattice constants, an edge state 
at the right-hand side of the sample with energy in the 
minigap decreases in energy, merges with the lower band, 
and later emerges as a left edge state. These calculations 
show that adiabatic ramping of the offset can produce a 
long-range, quantized, oscillatory mass current in a gen¬ 
eralized Fibonacci optical lattice. The effect does not 
depend on irrationality of the cut slope. This mass cur¬ 
rent could be detected, for example, by preferential load¬ 
ing of the edge states in a large-period lattice and direct 
imaging. Related effects have recently been observed in 
photonic waveguide lattices iHiiia, and recent theoretical 
and experimental work indicates that bulk atomic states 
can be pumped in a similar way in optical superlattice 
potentials [45M5] . The cold atom context, uniquely, en¬ 
ables realization of topological pumping in the presence 
of tunable interactions, and with variable adiabaticity. 
Such experiments would represent a controllable realiza¬ 
tion of Thouless pumping of edge states [35], and could 
provide a powerful tool for dynamical topological control 
of atomic wave functions. 

The availability of the “hidden dimension” quite natu¬ 
rally allows us to investigate another aspect of quasicrys¬ 
talline physics: the role of their soft modes. While of 
course the optical lattice potentials here are externally 
imposed and hence do not have true dynamical Gold- 
stone modes, we may simulate the effects of phonons and 
phasons by suitable manipulations of the lattice and cut- 
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FIG. 4. (Color online) Edge state topological pumping by 
phason driving. Top: Energy states near the first minigap 
versus offset of cutting beam along a lattice vector, at the 
Fibonacci cut slope. As the offset between the lattice and the 
cut beam is varied, left and right edge states cross the gap. 
Coloring of points indicates center of mass, with the same 
mapping as Fig.[^ Bottom: Variation of spatial probability 
amplitude of an initial edge state as offset is adiabatically 
varied. Position is normalized to 0 at the left edge and 1 at 
the right edge. 


ting lasers: for instance, shaking the lattice parallel or 
transverse to the cut corresponds to driving a phonon or 
a phason, respectively. Phasons have important but in¬ 
completely understood effects on thermal and electronic 
transport in real quasicrystals m- This is of interest not 
only for fundamental reasons, but also because of poten¬ 
tial technological applications of quasicrystals’ anoma¬ 
lous electrical and thermal transport characteristics. The 
influence of phasons is not understood in large part be¬ 
cause of the experimental difficulty of disentangling the 
effects of domain walls, crystalline impurities, and disor¬ 
der from those due to phason modes. A unique aspect 
of the generalized Fibonacci optical lattice is that it en¬ 
ables direct oscillatory driving of phason modes. Mea¬ 
suring the response of the system to driving such modes 
at variable frequency would constitute a new kind of lat¬ 
tice modulation spectroscopy, in which the modulation 
occurs in the higher-dimensional space from which the 
quasiperiodic lattice is projected. This capability, im¬ 
possible in other quasiperiodic systems, should allow un¬ 
precedentedly specific investigation of phason physics. 


V. CONCLUSIONS 


In conclusion, we have described a tunable quasiperi¬ 
odic optical lattice, presented calculations of the proper¬ 
ties of quantum gases in such a trap, and shown that this 
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generalized Fibonacci optical lattice will enable exper¬ 
imental realization of topological pumping and phason 
spectroscopy. Artificial quasicrystals such as those we 
propose allow the exploration of arbitrary quasiperiodic 
geometries, unrestricted by the laws of chemistry. Cre¬ 
ation of a fully tunable quantum quasicrystal would open 
the door to a large range of exciting experiments beyond 
those discussed here, including extensions of these tech¬ 
niques to higher-dimensional quasiperiodic lattices. The 
proposed realization of the cut-and-project construction 
allows tuning across the dimensional crossover from pe¬ 
riodic 2D lattices to quasiperiodic ID chains, enabling 
direct investigation of descendant quasiperiodic phases 
of well-studied correlated 2D systems. The unique tools 
of atomic physics can also enable new types of experi¬ 
ments: Feshbach tuning of the scattering length would 
allow exploration of the poorly understood role of in¬ 
teractions in quasicrystals [48], and time-varying poten¬ 
tials would enable dynamical experiments impossible in 
static lattices, such as phason spectroscopy. Experiments 
on quasiperiodic optical potentials may ultimately prove 
complementary to synthesis and characterization of solid 
and photonic quasicrystals, and could open another con¬ 
ceptual angle of attack on the problem of designing and 
predicting the properties of these complex materials. 
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Appendix: Effective one-dimensional model for 
harmonic lattices 

In this appendix, we describe how we may construct 
systematically more accurate approximations of the effec¬ 
tive one-dimensional potential for atomic motion parallel 
to the cut, order-by-order in perturbation theory. Our 
notation follows that in the main text. As a first step, 
we approximate the cutting beam potential Vc at leading 
order as a harmonic well, 

Vc{x,y) = -Ace ^(“o) 

« -Ac + ^rnuj]_r]_ 
where we have defined uj± = ^ 


We now perform a rotation of the coordinate system, 
{x,y) —)■ (r_L,r||). In the new coordinates, the Hamilto¬ 
nian is (ignoring an unimportant constant) 


H = 


2m 2m 


+ 


Vc{rj_) + UL(r_L,r||) 


where 


14(r±,r||) = 

- Al sin^ 

- Al sin^ 



(A.2) 


(A.3) 


is the 2D periodic lattice potential written in the ro¬ 
tated coordinate system. Note that we have allowed for 
anisotropic x- and y- coefficients; this will be relevant to 
the tunability of the lattice, discussed below. Assuming 
that Ai <C Ac, we Hrst solve for the exact eigenstates of 
the trap potential (the ‘subbands’, to borrow terminol¬ 
ogy from semiconductor physics); we will then compute 
the effects of the lattice potential on these eigenstates via 
perturbation theory. First, we write 


H = Ho + Hi 


(A.4) 


where 


and 


Ho = ^+ Vcir^) 
2m 

Hi = ^ + VLir^,rn). 


(A.5) 


(A.6) 


Hq has exact eigenstates '!/’n(^||U-L) of the form 

ipn{r\\,r±) = Xbee{r\\)(l>nir±), (A.7) 

where (jjnixx) is an exact eigenstate of the harmonic mo¬ 
tion along r±: 



+ Vc{r±) 


4 ‘n{r±) 


J_'), 


(A.8) 


with 



(A.9) 


and Xfree is any function. Note that we have deliberately 
chosen to include the kinetic energy along ry as part of 
the perturbation Hi, as this simplifies the calculations. 

Next, we add in the effects of the ry dispersion and 
the lattice potential Vl- From (A.31, we see that Vl 
mixes the transverse and longitudinal motion. It is con¬ 
venient to account for this by computing corrections to 
an effective potential for the longitudinal motion I4ff(ry) 
order-by-order in perturbation theory. In other words, 
we restrict the transverse motion to a specified subband 
(here, we take n = 0 for specificity, but similar argu¬ 
ments apply, mutatis mutandis, for any n) and compute 
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corrections due to virtual fluctuations to higher subbands 
order-by-order in perturbation theory. 

For the case of the n = 0 subband, by a straightforward 
calculation, this approach yields an effective Hamiltonian 
for motion along ry, 

Pfl 

^feff=^+Kfi(r||) (A. 10) 

with m* the effective mass and 


Kiff(r||) = 


_ 

n>0 


nhui± 


„ ( llUf (A.11) 


In the above expression, we have introduced a shorthand 
for the inter-subband matrix elements of the lattice po¬ 
tential, 

^L”(ni) = J dr± (j)* {rj.)VL{ri\,rj_)cl)r,{rj_). {A.12) 

It is evident that the perturbative expansion is controlled 
by the ratio ~ of the inter-subband matrix 

elements to the subband splitting. This parameter is 
straightforwardly tunable in the proposed optical real¬ 
ization of the generalized Fibonacci lattice. Note that 
the effective mass will be corrected at higher orders of 
perturbation theory, and that we may need to be careful 
about degeneracies introduced at higher orders; a com¬ 
plete analysis of the perturbation theory is beyond the 
scope of the present work. However, some key qualitative 
features of the perturbative series, modulo these possible 
complications, may be extracted simply by studying the 
form of the perturbation series, as we show in the next 
section. 


1. Fourier analysis of effective potential 



FIG. 5. (Color online) Fourier spectrum of effective ID poten¬ 
tial. The semi-logarithmic plot shows the relative amplitudes 
of Bragg peaks in units of recoil energy produced by perturba¬ 
tion theory up to second order. We have assumed = A^. 
Note that this is the analytically obtained Fourier transform 
of the potential itself, rather than the numerically obtained 
Fourier transform of the atomic density, which is shown in 
Fig. 3 of the main text. 


{±K ± K', ±2iF, ±2K'}, with the signs all chosen inde¬ 
pendently. The relative amplitudes between these peaks 
are trickier to evaluate analytically, but may be readily 
computed numerically as needed; however, from the fact 
that these peaks only emerge at second order, they are 
accompanied by a factor of ^ ^ /"‘,with x oi y 

chosen according to whether we obtain G from the lead¬ 
ing peaks by adding K oi K', respectively. See Fig. for 
a plot of the calculated Fourier spectrum. 

We see that we recover a complicated sequence of 
Bragg peaks as we go to higher orders in perturbation 
theory. At each order we will find Bragg peaks at higher 
harmonics, but these will be correspondingly at sup¬ 
pressed amplitude. Generalizing the line of reasoning 
above, we see that perturbation theory up to order JV 
yields a set of Bragg peaks 


We now turn to an analysis of Fourier components of 
the effective potential 14®. Consider the Vbo term. Since 
there is a single occurrence of Vl at this order, we see 
that its Fourier transform with respect to ry will contain 
only the harmonics present in Vl- It is easy to see that at 
this order, the Fourier transform has Bragg peaks only at 
G S {±K,±K'}, where K = ^ cos a and K' = ^ sin a 
(throughout, we ignore G = 0 peaks as they correspond 
to an unimportant uniform offset of the energy). The 
corresponding t4ff consists of a pair of harmonics whose 
minima are respectively at lattice spacings of A = 2 cos a 
and A' = ^^ . For tana = 2/{l + y/b) = the ratio 
of these spacings is indeed the golden ratio, A'/A = r. It 
is also straightforward to show that the relative ampli¬ 
tudes of the Bragg peaks is given by V±k oc , 

V±K' oc , where ^ = ^Jh/mujx is the char¬ 

acteristic ‘oscillator length’ of the trap potential. 

At second order, we have two occurrences of Vl in 
I4ff, leading to Bragg scattering at wave vectors G G 


G&{mK + nK'} (A.13) 


where TO,n are integers with \m\ |n| = A/", with ampli¬ 

tude 


/g 


(Ai)l-I {A^p 


exp 


1 

4 


{mK + nK)^ 


.(A.14) 


The above discussion should make it evident that, at 
least in principle, the Fourier transform of I4ff is char¬ 
acterized by a singular continuous spectrum, as long as 
the cut angle is such that tan a is irrational: in this case, 
there is no algebraic relation between iF, K', and thus the 
set {mK + nK', m,nG Z} has no smallest vector. How¬ 
ever, as computed above, and commented on in more 
detail below, many peaks for large |m|, \n\ will have ex¬ 
tremely small amplitudes and are for practical purposes 
absent. Nevertheless, this demonstrates that the gener¬ 
alized Fibonacci optical lattice is strictly more quasicrys¬ 
talline than the bichromatic potentials created to date. 
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and should be sufficient to explore a variety of relevant 
physical questions. 


2. Comparison to ‘classical’ cut-and-project 

The quick decay of Bragg peak amplitudes at higher 
orders is a consequence of the single-harmonic form of 
the lattice potential and the ‘soft’ projection imposed by 
the Gaussian beam. Traditional treatments of quasicrys¬ 
tals (see, e.g. [36]) differ from this in two ways: (i) they 
discuss a multiple-harmonic lattice, such as the one pro¬ 
duced by a periodic delta-function array; and (ii) they 
incorporate a singular cut (for instance, a step function). 
As the quasicrystalline potential is the product of these, 
the Fourier spectrum of the resulting ID quasicrystal is 
obtained by convolving the Fourier transform of the lat¬ 
tice, which consists of equal-amplitude peaks at recipro¬ 
cal lattice vectors, with the slowly decaying Fourier spec¬ 
trum of a step function (the slow decay is a consequence 
of the step-edge singularity). The resulting Bragg peaks 
have an amplitude that decays very slowly (as a power 
law) with their wave vector. Although this slow decay 
is known to play an important role in d = 3 quasicrys¬ 
tals [in], we note that our approach nevertheless retains 
the dominant features of quasiperiodicity, namely a set 
of incommensurate Bragg peaks. Indeed, even purely 
bichromatic lattices are known to have a localization- 
delocalization transition. 

As a Hnal comment, observe that the leading term in 
our perturbative approach is a direct transcription of cut- 
and-project: we could have obtained this by convolving 
the harmonic potential with the Gaussian that describes 
the atomic density transverse to the cut axis (this is 
the ‘tube’ that the harmonically confined atoms are re¬ 
stricted to.) However, the higher-order terms (that gen¬ 
erate higher Bragg peaks) are in fact non-classical: they 
emerge due to virtual fluctuations, and have no classical 
analogue. (Note that a corollary of the previous obser¬ 
vation is that the strength of higher-order Bragg peaks 
may be enhanced by the inclusion of higher harmonics 
in Vl, as these would then contribute already at leading 
order.) 


3. Relation to Aubry-Andre potential 


A key feature of our set-up is the inherent tunabil- 
ity afforded by the cut-and-project approach. As an ex¬ 


ample, we sketch a prescription of how to achieve the 
Aubry-Andre limit in the cut-and-project approach. To 
do this, we exploit the ability to impose an anisotropic 
lattice potential in two dimensions. In the limit = 0, 
we simply have a stripe modulation in 2D. The effective 
ID potential is periodic, with spacing given hy Xl/ cos a. 
For a strong Gaussian beam Ac ^ we can restrict 
ourselves to leading order in perturbation theory, so that 
it suffices to consider the leading harmonic of this ID 
potential, with strength Vb oc V±k oc A| (here and be¬ 
low, we have ignored factors of 0 ( 1 ) as we are interested 
in the scaling rather than precise numerical factors). As¬ 
suming a relatively deep lattice, we may approximate the 
projected ID lattice by a tight-binding chain: in second- 
quantized form, we have 

Fleff « (A. 15) 


with a hopping matrix element (using standard tech¬ 
niques, m) 


j 





(A.16) 


where En is the effective ID recoil energy (this is pro¬ 
portional to the recoil energy of the 2D lattice, but as 
is the case for other quantities, may differ by factors of 
0 ( 1 ) due to the projection.) 

Next, we impose a weak periodic potential in the y 
direction, by allowing 7 ^ 0. If we assume that 
A)) A\ (A)))^/?iaj_L, we may simply consider the 
additional contribution to the effective potential due to 
as a weak perturbation to the dominant contribu¬ 
tion due to A)( (while still ignoring transitions to higher 
subbands of the trapping beam). This additional con¬ 
tribution comes from harmonics at K' = AT tan a, with 
strength Vi oc V±k' oc Since A^ :S> A^, this may 
be incorporated as a weak onsite modulation to the ID 
tight-binding Hamiltonian, 


HeS ■ 


i 

+ + U cos{K'ri))blbi 


(A.I7) 


which is equivalent to that of the Aubry-Andre model 
with a phase offset of zero. We may recover the full 
Aubry-Andre mode l by also including an offset phase in 
the A^ term of Eq.(A.3). 
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